ANALYSIS AND SIMULATION OF 
OF WIND-DRIVEN CURRENTS DURING 
THE MIXED LAYER EXPERIMENT 
(MILE) 



Jose M. Fernadez Lopez 



NAVAL POSTGRADUATE SCHOOL 

Monterey, California 




THESIS 



ANALYSIS AND SIMULATION OF WIND-DRIVEN CURRENTS 
DURING THE MIXED LAYER EXPERIMENT (MILE) 


by 




Jose M. Fernandez Lopez 




March 1981 




Thesis Advisor: 


R. W. Garwood 



Approved for public release; distribution unlimited 



SECuniTV CLAISiriCATlON OP THIS PAGE (Whm% Dmf Bntmfmd) 



REPORT DOCUMENTATION PAGE 


READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


1. PEPOPT NUMPCP 


2. GOVT ACCESSION NO. 


3. PECIPIENT’S CATALOG NUMaCR 


4 . title (M%d Subtttim) 

Analysis and Simulation of Wind-Driven Currents 
During the Mixed Layer Experiment (MILE) 


S. TYPE OF REPOPT 4 PERIOD COVERED 

Master's Thesis 
March 1981 


4. performing OPG. PEPOPT NUMSEP 


7. AuTMOPr#; 

Jose M. Fernandez Lopez 


4 . contract op grant NLMOCRr*; 


f. PEPfOPMING OPGANIZATION NAME AND AOOPCSS 

Naval Postgraduate School 
Monterey, California 93940 


10. PROGRAM element. PROJECT TASK 
AREA 4 WORK UNIT NUMBERS 


n. CONTMOLLINC OfftCE NAMC AND AOOMCSS 

Naval Postgraduate School 
Monterey, California 93940 


• 2. MEPOnr DATE 

March 1981 


11. number oe pages 

52 


u. MONITOPING agency NAMC 4 AO 0 m €%%( H dlUbrwnt irom Cot%tr«tting OHIcm) 


IS. security class. (oI ihl» r.Rero 

Unclassified 


Mo. DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


16. DISTPiauriON STATEMENT ro/ P«porO 

Approved for public release; distribution unlimited. 


17. OISTPI0UTION STATEMENT (ot thm mtmtfmci in Block 20» U dlllmtofit from Boporl) 


It. supplementapy notes 


19. KEY WOPOS rConflnu* on toworcm oldo II noeoomofr ond tdonilly ky block numborj 

Wind-Driven Current 

Mixed Layer Experiment (MILE) 

One-dimensional Current Model 


20. ASST PACT cConllnuo on rovorom oldo II nmcomoocy md Idontltf by block mmiborj 

The wind stress calculated from wind velocities measured during the Mixed 
Layer Experiment (MILE) was used as input to a one-dimensional wind-drive 
current model. These model results are compared with observed currents 
from the MILE-1 bouy, showing a qualitative agreement. 

MILE was an examination of the upper ocean carried out near Ocean Weather 
Station Papa during a 20-day period, August and September of 1977, which 
was characterized by two major wind events. 



DO 1473 EDITION OP I NOV tt IS oatOLETE 

S/N 0102-014-««0l I 



1 



iCCUfllTY CUAfSiPiCATlON OP THIS PAGE (Wht% Ddf 




f 








cu 















The observed currents have been analyzed to obtain information about 
their behaviour that could be used in the tuning of the model. For a 
simulation of the entire period the results are considered only 
satisfactory. 



i 



3 Fornj 1473 2 

0*^2-014-6601 



iicuatTv ckAMirtCATioM 0^ ot» #<*••»•*» 



Approved for public release; distribution unlimited 



Analysis and Simulation of Wind-Driven Currents 
During the Mixed Layer Experiment (MILE) 



by 



Jose M. Fernandez Lopez 
Lieutenant Commander, Spanish Navy 



Submitted in partial fullfillment of the 
requirements for the degree of 

MASTER OF SCIENCE IN OCEANOGRAPHY 

from the 

NAVAL POSTGRADUATE SCHOOL 
March 1981 



ABSTRACT 



The wind stress calculated from wind velocities measured during 
the Mixed Layer Experiment (MILE) was used as input to a one-dimensional 
wind-driven current model. These model results are compared with ob- 
served currents from the MILE-1 bouy, showing a qualitative agreement. 

MILE was an examination of the upper ocean carried out near Ocean 
Weather Station Papa during a 20-day period, August and September of 
1977, which was characterized by two major wind events. 

The observed currents have been analyzed to obtain information about 
their behaviour that could be used in the tuning of the model. For a 
simulation of the entire period the results are considered only 
satisfactory. 
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I. INTRODUCTION 



A. PURPOSE OF THE STUDY 

The objective of this research was to investigate how well an one- 
dimensional model of the wind driven current could explain the observed 
near surface current during the Mixed Layer Experiment (MILE). Measure- 
ments of the near-surface currents were taken during MILE, August- 
September 1977, in the vicinity of Ocean Weather Station (OWS) PAPA. 

During this period, the data show the response of the upper ocean to the 
passage of two wind events. 

This study investigates the hypothesis that surface wind can account 
for a large part of the energy and variability of inertial oscillations 
near the ocean surface and that these oscillations are predominantly 
locally generated, as was demonstrated by Pollard and Millard (1970) at 
another site. 

To accomplish this task, a one-dimensional wind-driven current model 
was solved numerically using a modified "leapfrog" scheme. The modifi- 
cation involved employed current values, u and v, which were vertically 
averaged at each time step to a depth equal to that of the observed mixed 
layer, simulating strong boundary layer mixing. The mixing depth used 
was a time-average of the mixed layer depth obtained from buoy temperature 
records and verified with Plessey CTD data taken from the R/V 
OCEANOGRAPHER. 

The data used to check the response of the model were those obtained 
by nine near-surface current meters, at depths from 5 to 32m, of the 
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MILE-1 Mooring. The wind stress was computed from wind speed and di- 
rection as recorded by the R/V OCEANOGRAPHER during the experiment. 

It was hoped that the non-wind-driven current would be sufficiently 
small and invariant so that the hypothesis could be tested. 

In the next section, the current meter observations are analyzed by 
using a Fast Fourier Transform. The power spectra of currents from three 
significant depths illustrate how the currents are dependent on the 
surface wind stress. 

In section three the equations on which the model was based are given. 
It is shown how from those equations a numerical solution is obtained to 
solve the initial value problem that we are dealing with and how the 
model results were computed. 

In the last section the results and conclusions are given. 

If there is no other indication, the abscissa in all figures gives 
time in hours from the start of the experiment, 0500 GMT hours August 
18th. There are 460 hours of recorded currents in all. 

B. IMPORTANCE OF THE PROBLEM 

As Csanady [1981] explains: "The problem of circulation is to de- 

scribe and understand the pattern of the residual or longer-term water 
particle displacement. The distribution of important water properties, 
e.g., temperature, salinity (...), depends critically on the pattern of 
circulation." The capability of modeling wind-driven currents could 
help to explain one component of that pattern of circulation. 

Modeling wind-driven current is needed in order to improve upper ocean 

thermal structure models, which have great importance in ASH. Fisheries, 

search and rescue operations, and control of pollution in the ocean are 

other fields where wind-driven currents must be considered. 
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C. LITERATURE REVIEW 



Horizontal inertial currents occur in thin layers with a relatively 
reduced horizontal extent, Webster [1968] concludes that inertial 
oscillations are essentially transient phenomena. For a fixed location 
inertial oscillations at different depths may have diverse origins. 
According to Pollard and Millard [1970], Gonella [1971] and Kundu [1976], 
for the surface layer of the ocean, the wind stress is the principal 
mechanism for their generation. 

Inertial currents rotate clockwise (Northern Hemisphere) with a 
period T = 7r/{lfij sin 6), where is the earth's rotation vector and 
9 is the latitude of the observation point. The transmission of the 
movement from the surface to deeper layers may normally be carried out 
either by turbulence in an unbounded homogeneous medium, or by boundary 
effects in a strongly stratified medium. 

To test the dependence of inertial oscillations as a function of 
wind stress the theory of Ekman [1905] was applied in an impulsive 
system. In this theory, the wind provides a stress on the sea surface, 
and the motion is viscously transmitted to lower layers. 

Figures 1 and 2 show four typical examples of the measured vertical 
structures. In figure 1 both temperature and velocity profiles are con- 
tinuously stratified: the negative gradients of temperature show that 

there is no mixing. However, the profiles in figure 2 show a change in 
temperature of less than 0.2° C in the upper 28m, indicating that the 
water is nearly homogeneous and the current is also nearly uniform. The 
more the ocean is stratified, the greater is the tendency for agreement 
with Ekman 's model which was developed for a homogeneous ocean. As 
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Figure 1. Stratified profile. 
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Gonella [1971] concludes from his work in the Mediterranean this agree- 
ment comes about "because the assumption of a constant coefficient of 
eddy viscosity is more realistic in a continuously stratified ocean than 
in an entirely homogeneous ocean where the turbulence may be too intense". 

The approach to the problem in this study is similar to the "slab" 
model of Pollard and Millard [1970] and Garwood [1976], using a depth 
integral of the momentum equation to solve the problem of turbulent 
momentum transfer within the boundary layer. 

The essential principle is that the integrated inertial force, in- 
cluding Coriolis and local time derivative accelerations, just balances 
the applied surface stress less a small damping term. 

To interpret actual current records, two modifications were intro- 
duced in the present work: (i) the momentum equation was integrated to 

a depth equal to that of the mixed layer in each time step, and (ii) a 
linear damping term was applied. By the former procedure, the modeled 
mixed layer becomes unstratified with f = 2 tt/T as its only natural fre- 
quency. The latter step has the effect of changing that frequency to a 
slightly smaller value. This may not be a good property when modeling 
the ocean, where inertio-gravity waves have frequencies betv;een f and 
N, the Brunt-Vaisala frequency. Since N is usually greater than f, 
the observed near-inertial currents frequently have a frequency slightly 
greater than f, [Pollard and Millard, 1970]. 

In the one-dimensional model used here surface winds are expected to 
account for a large part of the variability of the energy of the near- 
inertial oscillations near the surface. It was also expected that, as 
in other models [Pollard and Millard, 1970; Kundu, 1976], the changes in 
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phase and amplitude of the forced oscillation could be simulated with 
some success. 

Another problem to be considered is how the energy at inertial fre- 
quencies is removed from the near-surface layers. Pollard [1970], 

Kroll [1975], and Kase [1979] show that one possibility is that energy 
is dispersed through the seasonal thermocline by transferring into in- 
ternal waves in the near-inertial band. Other removal mechanism could 
be the wind [Pollard and Millard, 1970], and in the next section 
evidence will be presented that shows how the wind by, shifting its 
direction, may remove part of the kinetic energy that had been accumu- 
lated earlier. 

Pollard [1980] shows that the amplitude of inertial oscillations can 
fall to zero within a few tens of kilometers, and that the decay of energy 
is due to a slow downward dispersion associated with horizontal spreading, 
suplemented by the effect of changing winds. In the present work there 
is particular interest in the downward dispersion and the effect of 
changing wind direction. 
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II. THE MILE DATA 



A. DESCRIPTION 

From 19 August to 6 September 1977, two surface moorings were de- 
ployed in the vecinity of OWS PAPA. The MILE-1 mooring, which provided 
the data used in this study was located at 49° 37 'N and 135° 6'W. 

This mooring consisted of an 8m surface toroid anchored with 0.9 scope 
in 4360m of water. Shackled in the mooring line were 19 VACM's of 
which the upper nine were used in this study; those at mean depths of 
5, 8, 11 , 14, 20, 23, 26, 29 and 32m. 

The reason for choosing only those current meters was that the 
mixed layer was at all times more shallow than about 30m, so the measure- 
ments below 32 meter were not considered to be of much significance in 
testing the hypothesis. 

All recorded data were pre-processed and provided on magnetic tape 
in intervals of 15 minutes for a 19-day period, with the following 
exceptions ; 

At the 5m level , there are only 5 days of current record. 

At the 8m levle, there are only 16 days of current record. 

At the 14m level there are only 3 days of current record. 

Temperature records at all levels are complete. 

Wind data were recorded hourly aboard the R/V OCEANOGRAPHER. 

Figures 3 and 4 show those wind values. Two strong wind events can be 
seen one strong storm in 22-24 August (starting time 100 hours), and a 
weaker storm in the 31 August-1 September period (starting time 310 
hours). The six-hourly sea-surface pressure maps reveal that both storms 
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Figure 3. Wind speed vs time. 
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Figure 4. Wind speed vs time, u and v component. 
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were low pressure disturbances which developed rapidly in the Gulf of 
Alaska and moved toward land in an East-Northeast direction. The center 
of the low on 31 August passed directly over station P, causing calm 
conditions for a few hours. That is shown by the double pulse in figure 3. 

B. ANALYSIS 

Figure 5 shows hourly temperatures at the studied levels and illus- 
trates the variability of the temperature structure of the layer above the 
seasonal thermocline. Homogeneous layers are indicated in this figure when 
temperature traces from two or more depths coincide. During stratified 
conditions, the vertical or horizontal movement of the water column past 
the sensors causes apparent temperature variability. Most of this vari- 
ability appears to be due to tidal -period internal waves. Above the 
30 meters level, homogeneous and stratified conditions alternate in re- 
sponse to atmospheric conditions, as can be seen in the power spectra in 
figure 5a. 

With the rapid increase of the wind early on 22 August, at about 95 
hours on the time axis, the temperature at 5m dropped and within four 
hours the upper 23m of the water column became vertically homogeneous 
due to strong mixing. During this strong wind event, the mixed layer 
deepened to nearly 35m. After the storm died, by 24 August stratifi- 
cation was reestablished between 17m and the base of the earlier 30m mixed 
layer. During the second storm the weaker wind event of 27-28 August 
mixed the upper layers to about 22m (figure 6). 

Figures 7, 8 and 9 are the power spectra for the u and v com- 
ponents of the observed current at three different levels, 11, 23 and 
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Figure 5a. Power spectrum of the temperature at two levels. 
Frequency: c.p.h. 
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Figure 6. Mixed layer depth variation with time, as computed 
from MILE-1 temperature data buoy. 
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32m, for three different time windows. In figure 7 the time window is 
from 0 to 260 hours, only the effects of the first storm are present. 

In figure 9, the time window is from 200 to 460 hours, and only the 
second storm is present. For the spectrum of figure 8, the time window 
is from 100 to 360 hours, and the time period includes both storms. 

Comparison between different depth levels in figures 7 and 9 indi- 
cates different amounts of kinetic energy for each depth at near-inertial 
frequencies. The first storm, as can be seen in figure 4, was charac- 
terized by a pronounced shift from steady easterly winds to strong 
westerly winds. That strong change after a period of at least 100 hours 
of steady winds could be the means by which part of the energy transfer 
between the ocean and the atmosphere could have eliminated the existing 
inertial current before generation of a new inertial current. 

Since the upper ocean is assumed to be one-dimensional there must 
be some unsteady mechanism by which the vertical transfer of momentum 
is modulated in time, causing the observed phase and magnitude differences 
between the 32 meter level current and the surface. 

The response of the current in the upper layers to the second storm 
was quite different from the response to the first storm. In this latter 
case the wind shifted from having a northerly component to a strong 
southerly one and, in a few hours, to a northerly component again 
(figure 4). The upper level currents initially appear to be attenuated by 
the southern component winds but then are reinforced when the wind di- 
rection shifts toward the south. Also these winds have a smaller inten- 
sity and duration, and as it is shown by figure 6 the mixed layer did not 
deepen as much as in the first storm. 
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Figure 7. Power spectrum of observed current. 

Time window 0 to 260 hours. 

(Abscissa: c.p.h.; ordinate: (cm/s)'^). 

1 : inertial frequency. 
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Figure 8. Power spectrum of observed current. 

Time window : 100 to 360 hours. p 
(Abscissa: c.p.h.; ordinate: (cm/s)^). 

1 : inertial frequency. 
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Figure 9. Power spectrum of observed current. 

Time window : 200 to 460 hours. ^ 
(Abscissa: c.p.h.; ordinate: (cm/s) ). 

1 : inertial frequency. 
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Because the mixed layer had undergone a shallowing process as dis- 
cussed earlier, the transfer of energy to the 32m level was hindered by 
stratification. 

Figures 10 and 11 are the power spectra for the u and v compon- 
ents of the current for the total period of observation. The spectra 
for 32m level exhibit a double-peaked feature at the near-inertial period. 
As Pollard [1970] suggests, at that depth it is unlikely that the ampli- 
tude of inertial oscillations can be generated by the action of a single 
storm. 

The time series of observed currents are shown in figures 12 to 19. 

In those figures, the plotted running mean, an average over one inertial 
period, shows the non-inertial current at each level. It should be 
recognized here, that an undetermined part of this non-inertial current 
may be quasigeostrophic and not directly related to the local wind 
conditions. 
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Figure 10. Power spectrum of observed current. 

Total period (u-component) . ^ 

(Abscissa: c.p.h.; ordinate: (cm/s)'^). 

1 : inertial frequency. 
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Figure 11. Power spectrum of observeti current. 

Total perio(i (v-component) . ^ 

(Abscissa: c.p.h.; or(Jinate: (cm/s)^). 
1 : inertial frequency. 
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Figure 12, Observed current, total period Cu-component) . 

Running mean (*), average over 1 inert, period. 



28 



CM/S (V-COMP.) CM/S (V-COMP.) CM/5 (V-COMP 



43 . 




- 40 . I ! ! I : — J i 1 1 I ^ 

0 . 100 . 200 . 300 . 400 . 500 . 

HOURS, (DEPTH 11 M.) 




HOURS. (DEPTH 23 M.) 




HOURS. (DEPTH 32 M.) 



Figure 13, 



Observed current, total period (v-component) . 
Running mean {*), average over 1 inert, period. 
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Figure 14. Observed current: 0 to 260 hours (u-component) . 

Running mean (*), average over 1 inert, period. 
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Figure 15. Observed current: 0 to 260 hours ( v-component) . 

Running mean (*), average over 1 inert, period. 
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Figure 16. Observed current: 100 to 360 hrs. (u-component) . 

Running mean (*), average over 1 inert, period. 
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Figure 17. Observed current: 100 to 360 hrs. (v-component) . 

Running mean (*) , average over 1 inert, period. 
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Figure 18. Observed current; 200 to 460 hrs. (u-component) . 

Running mean (*), average over 1 inert, period. 
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Figure 19. Observed current: 200 to 460 hrs. (v-component) . 

Running mean (*), average over 1 inert, period. 
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III. DESCRIPTION OF THE MODEL 



A. THEORETICAL EQUATIONS 

In this study, one-dimensional conservation of momentum in the mixed 
layer is assumed. This conservation of momentum and the condition of 
incompressibil ity are reflected by the Navier-Stokes equations of motion. 







Vp + ^ + vV 



2 ^ 

V 



8u. 

1 



8x, 



= 0 



setting: u = u + u + u' 

g 

V = V + V + v' 

g 

where u and u' are the horizontal mean and fluctuating components of 
x-component of non-geostrophic velocity, and u^ is its geostrophic 
component; similarly v, v' and v^ are the corresponding y-components. 
With the Boussineq approximation and Reynolds averaging, 



(Ug + 0 u') = f(Vg + 5 + V) - 



3u'w' 3u' v‘ u' u' 

3z ~ 3y ~ 3x 



+ V 



I 2 2 2 1 

3 u , 3^v , 3“^w 
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(Vg + ; + V) = -f(Ug t s + u') - 



3v'w' _ 3v ‘ v' _ 3u'v‘ 
9z " 3y ~ 3x 



+ V 



12 2 2 1 

3 u , 3 V , 3 w 

3x‘^ 3y^ 3z J 



Assuming 

(a) Reynolds stresses are much greater than viscous stresses. 

(b) Vertical velocity field vanishes, w=0. 

(c) Geostrophic flow is separable from wind driven flow, 
and neglecting the horizontal density gradient 



fu 



g 



I_ 12. 

Pq 3y 



fv 



g 



_ 1 3p 
- P(, 3x 



Then the equations for wind-driven flow become 



li = 

3t 



fv - T— u'w’ 

oZ 



[ 1 ] 



3v 

3t 



fu - ^ Vw' 



[ 2 ] 



B. NUMERICAL METHOD OF SOLUTION 

To solve equations [1] and [2] a vertical grid from the surface to 
35m, was used. The values of u(z,t) and v(z,t) were then solved 
numerically as an initial -value problem. 
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Taking the equations for wind-driven flow and applying an eddy 



viscosity closure. 



-u'w' = kf 



— I — r f 3v 

-Vw' = k ^ 

and assuming k constant, gives: 



f = f5 . k g [3] 



8V ^ 



The boundary conditions are: 
a. At the surface: 



[4] 






’0 



3z 



-^''^0 = 3z 



3 ^ 0 ) _ 



b. 



where t (t) and x (t) are wind 
X y 

directions respectively at time t 
of sea water. 

At z = -D < -h, a "slip" condition 



u'w'(-D) = 0 



v'w'(-D) = 0 



stress in x and y 
and p is the density 

is prescribed. 
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With the eddy viscosity closure. 



3u 

3z 



(-D) = ||- (-D) = 0 



The numerical model used to simulate the current was based on a 
"leapfrog" scheme, with m indicating space (depth) and n time, then 
u(z,t) = u(mAz,nAt) = u(m,n). In this manner the component of equations 
[3] and [4] can be written: 



3u _ u(m,n+l ) - u(m,n-l ) 
3t 2At 



3v _ v(m,n+l ) - v(m,n-l ) 

3t ~ 2At 

fv = fv(m,n) 
fu = fu(m,n) 

3^11 k 

k = y (u(m-l ,n-l ) - 2u(m,n-l ) + u(m+l ,n-l ) ) 

3z^ (Az)^ 

3^v k 

k = o (v(m-l ,n-l ) - 2v(m,n-l ) + v(m+l ,n-l ) ) 

3z^ (Az)^ 

Solving for the (n+1) values, 

u(m,n+l )=Ru(m-l ,n-l )+(l-2R)u(m,n-l )+Ru(m+l ,n-l )+2fAtv(m,n) 
v(m,n+l )=Rv(m-l ,n-l )+(l -2R) v(m,n-l )+Rv(m+l ,n-l )-2fAtu(m,n) 
where: R = 2kAt/(Az)^. 

At the surface m = M and application of the boundary conditions 
gives : 
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u(M,n+l) = 2Ru(K-1,n-l) + ( l-2R)u(M,n-l ) 



+ 2fAtv(M,n) + 2RAz Tx(n)/(kp) 



v(M,n+l) = 2Rv(M-l,n-l) + (l-2R)v(M,n-l ) 
- 2fAtu(n,n) + 2RAzTy(n)/(kp) 



The numerical model is stable when the value of R is equal to or 

less than 0.5. Assuming an eddy viscosity coefficient k of 23 cm/s 

and an increment of depth of Az = Im, and a time step At = 100 

seconds gives an R equal to 0.46. 

The initial conditions were established as those observed values of 

the current speed for each component, u and v, at time 0500 GKT 

August 19th, at the first recorded peak of the v-component of the current. 

The surface stresses tx and ly were computed every hour as a 

function of the hourly wind speed data for each component w and w , 

X y 

usi ng 



2 c 

T = W (W^ + W^) 
X x' X y' 



1/2 






1/2 



T = W (W^ + W^) Cr,P 

y y X 'j' D^a 



,-3 



with a value for the drag coefficient of 1.3 x 10 

Values for the u and v components were calculated for each 
level, from the surface to 35m at intervals of Im, at time steps of 
100 seconds. 

Every hour a new input of wind stress was applied, setting as 
initial conditions for the next run an integrated value for depth in the 
mixed layer. 
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Mixed layer depth as a function of time was deduced from the Mile-1 
temperature time series and checked with Plessey CTD data taken from the 
R/V OCEANOGRAPHER. 

The model was initially run without any damping term, and the 
results were acceptable until the onset of the second storm. At that 
time a clear shift in phase between the recorded and modeled currents 
was observed. In an attempt to avoid a steadily worsening phase problem 
a linear damping term was introduced. According to Pollard and Millard 
[1970], this decay factor models the dispersion effect. Then the 
equations for wind driven flow become: 

II = rt . |_ (JV ) - cG 

8v r'* ^ / i T\ “ 

= -fu - gj (Vw') - cv 

-5 -1 

A damping coefficient of c = 5.03 x 10 sec was chosen from 
Pollard [1970]. Thsi gives an e-folding decay time of 2.3 days. 
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IV. RESULTS 



Figure 20 through 22 show the modeled currents compared with the 
recorded ones at three diffe rents levels. 

As was discussed earlier, the dominant frequencies of the model 
and the ocean are slightly different, see figures 23 and 24. There- 
fore when the wind decreases and the current is almost totally inertial 
under small forcing the modeled and observed oscillations will drift 
out of phase. This is most apparent in the period before the second 
storm, when the wind speed dropped to a minimum. 

The model does not have the capability to reproduce the observed 
current driven by the sharp change in directions of the wind during the 
second storm. The modeled and observed currents become out of phase 
at about time 330 hours. After the wind again becomes more steady in 
direction, the modeled and observed currents again agree in phase. 

Other discrepancies between observed and modeled current could be 
caused either by observational errors or by failure of the model to 
properly treat the physics of the problem, such as neglecting advection. 

It is clear that the dissipation of mean kinetic energy is too much 
in the model using R = 1/2.3 days. Figure 22 shows that at the 32m 
level the modeled current is stronger than the recorded one. 

At the other levels the opposite happens. The model needs to be 
tuned to get a better amplitude comparison in the mixed layer. 

Doing this tuning by only comparison with observed currents is 
difficult because of the probable errors in the measurements. As 
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Figure 20. Comparison between observed and modeled currents. 
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Figure 21. Comparison between observed and modeled currents. 
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Figure 22. Comparison between observed and modeled currents. 
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Figure 23. 



Power spectrum of modeled current (u-component) . 

2 

(Abscissa: c.p.h.; ordinate: (cm/s) ). 
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Figure 24. Power spectrum of modeled current ( v-component) . 
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suggested by Halpern et al 1^981], there is spurious rotation of the 
rotor of the current meters produced by motions of the surface- 
following buoy mooring. 

Tuning could be done with a different treatment of the eddy vis- 
cosity as v/ell as alternate parameterization of the mean kinetic energy. 
The first approach would be possible when turbulent viscosity of the 
water as a function of depth and stratification is better known, per- 
mitting detailed quantitative calculations of the vertical stress 
profi le . 

Two other factors that could have a relative important effect on the 
currents and are not included due to the limited scope of this study is 
the process of inertial-gravitational wave propagation and the resultant 
horizontal dispersion of mean kinetic energy. 
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